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PREFACE 

A  difficult  and  recurrent  problem  in  many  applications  of  fluid 
mechanics  is  the  description  of  the  convection  of  a  passive  scalar 
quantity  in  a  turbulent  flow  field.  Examples  of  the  diverse  areas  in 
which  this  problem  arises  are  in  reentry  physics,  when  electron  con¬ 
centration  is  effected  by  wake  turbulence,  or  in  pollution  studies, 
when  the  passive  scalar,  pollutant  concentration,  is  transported  by 
the  turbulent  flow  field  of  the  surrounding  medium. 

A  theoretical  investigation  of  this  problem  is  presented  in  this 
Memorandum,  using  a  systematic  approximation  based  on  an  expansion  in 
Weiner-Hermite  functionals.  This  approximation  leads  to  analytic  so¬ 
lutions  to  two  crucial  problems  in  the  theory  of  turbulent  diffusion: 
the  diffusion  of  a  passive  scalar  from  a  point  source,  and  the  spec¬ 
trum  of  a  statistically  homogeneous  scalar  --  both  in  a  homogeneous 
turbulent  flow. 

The  research  described  here  is  part  of  RAND’s  work  on  reentry 
aerodynamics  for  the  Advanced  Research  Projects  Agency. 

The  author,  P,  G.  Saffman,  is  a  Professor  of  Fluid  Mechanics 
at  the  California  Institute  of  Technology,  and  a  consultant  to  The 
RAND  Corporation. 
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SUMMARY 


The  velocity  and  concentration  fields  in  a  homogeneous  turbulent 
flow  convecting  a  passive  scalar  are  expanded  in  series  of  Wiener- 
Hermite  functionals.  The  equations  of  motion  lead  to  integro-dif f er- 
ential  equations  for  the  coefficients  of  the  expansion.  The  expansion 
and  the  equations  are  truncated  at  the  first  stage  that  gives  nontriv¬ 
ial  results,  and  analytical  solutions  are  obtained  for  the  concentra¬ 
tion  field  with  the  velocity  field  given  and  with  zero  molecular  con¬ 
ductivity. 

The  diffusion  from  a  point  source  is  considered  first.  It  is 
shown  that  the  mean  concentration  satisfies  an  integro-dif ferential 
equation  with  diffusive  and  wave- like  properties,  which  reduces  for 
large  time-values  to  the  diffusion  equation  in  accordance  with  the 
classical  theory  of  G.  I.  Taylor.  The  effective  diffusivity  is  deter¬ 
mined  in  terms  of  Eulerian  velocity  correlations. 

The  spectrum  of  a  statistically  homogeneous  distribution  of  the 
scalar  is  considered  second.  It  is  shown  that  if  the  spectrum  is  ini¬ 
tially  concentrated  around  the  wave  number  of  the  energy-containing 
eddies,  an  equilibrium  is  set  up  for  larger  wave  numbers,  with 

a  2 

F (k)  =  -2-  E (k) 

3u 

where  f(k)  and  E(k)  are  the  concentration  and  energy  spectra,  respec- 
“2  2 

tively,  and  0  and  u  are  the  mean-square  concentration  and  the  mean- 
square  velocity  component.  Combined  with  a  modification,  to  take  ac¬ 
count  of  small-scale  intermittency ,  of  Batchelor’s  theory  of  the  spec¬ 
trum  for  very  large  wave  numbers,  this  equation  is  used  to  predict  the 
decay  rate  of  the  concentration  intensity. 
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I.  INTRODUCTION 

The  diffusion  and  mixing  of  convected  quantities  in  a  turbulent 
flow  is  a  subject  of  great  importance  and  has  attracted  much  attention. 
For  a  passive  scalar  in  an  incompressible  fluid  with  uniform  proper¬ 
ties,  the  concentration  9(x,  t)  satisfies  the  equation 

+  u  •  V0  =  kV^Q  (1) 

at  — 

where  u(x,  t)  is  the  velocity  field  and  k  is  the  molecular  diffusivity. 
Given  the  velocity  field  u,  Eq.  (1)  is  linear,  unlike  the  nonlinear 
Navier-Stokes  equations  for  the  velocity  u  itself.  But  this  advantage 
is  apparent  rather  than  real,  and  a  linear  equation  with  random  coef¬ 
ficients  seems  to  present  as  hard  a  problem  as  a  nonlinear  one. 

The  problem  has  two  distinct  aspects.  First,  there  is  the  ques¬ 
tion  of  how  an  inhomogeneous  distribution  is  mixed  and  transported  by 
the  fluid  motion.  This  can  be  subdivided  further  into  diffusion  from 
a  point  source  (sometimes  called  one-particle  analysis),  first  studied 

by  Taylor, ^  and  the  spreading  of  a  cloud  (two-particle  analysis), 

(2) 

first  investigated  by  Richardson. v  A  survey  of  this  work  is  to  be 

(3) 

found  in  an  article  by  Batchelor  and  Townsend.  There  is  also  the 

difficult  matter  of  deciding  to  what  extent  the  mixing  due  to  turbu¬ 
lence  can  be  described  by  an  eddy  diffusivity,  and  how  the  eddy  diffu¬ 
sivity  (if  it  is  a  reasonable  approximation)  is  related  to  the  statis¬ 
tics  of  the  velocity  field.  The  second  aspect  concerns  the  intensity 
decay  and  spectrum  of  a  homogeneous  distribution  in  an  approximately 
homogeneous  flow.  The  pioneering  work  on  this  subject  was  done  by 

Obukhov  and  Corrsin,  and  a  unified  treatment  has  been  given  by 

(4) 

Batchelor . v 

Although  intensive  effort  on  these  problems  has  uncovered  many 
details  of  the  processes  and  provided  not  too  unsatisfactory  empirical 
laws  for  practical  applications  in  certain  cases,  a  basic  understanding 
of  the  fundamental  questions  is  still  seriously  incomplete.  Part  of 
the  difficulty,  but  by  no  means  all,  lies  in  ignorance  about  the  veloc¬ 
ity  field  itself. 
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For  diffusion  from  a  point  source  in  an  isotropic  flow,  Taylor 
showed  in  effect  that  for  large  values  of  time  the  mean  concentration 
is  described  by  the  heat  diffusion  equation 

||  +  u  •  ve  =  h  v2e  (2) 

ot  —  e 

where  U  is  the  constant  mean  velocity  and  is  an  effective  diffusiv- 
ity,  related  to  the  Lagrangian  auto-correlation  of  the  velocity  follow¬ 
ing  a  fluid  particle.  The  definition  and  relation  are 

O  00 

U  VT)  =  (u(t)u(t  +  t)^  J  Ke  ~  u2  [  Rl(t)  <3t  (3) 

where  u(t)  is  the  component  of  velocity  fluctuation  of  a  fluid  particle 
and  u  is  the  root  mean  square.  One  of  the  unsolved  problems  of  turbu¬ 
lence  is  the  relation  between  the  Lagrangian  auto-correlation  and  the 
more  readily  measured  Eulerian  correlations,  although  crude  theoreti¬ 
cal  estimates  have  been  made  in  prof usion. ^ 9 ^ 

All  previous  work  on  the  spectrum  and  decay  of  a  homogeneous  dis¬ 
tribution  has  been  based  on  the  hypothesis  of  a  cascade  and  the  in¬ 
dependence  of  large  and  small  wave  numbers,  and  has  been  closely  re¬ 
lated  to  the  Kolmogorov  universal-equilibrium  theory.  However,  the 
Kolmogorov  theory  itself  has  been  both  modified^)  and  questioned^^ 
in  recent  years,  and  it  seems  that  the  evidence  in  its  favor  is  by  no 
means  as  strong  or  as  convincing  as  was  once  supposed.  Thus  it  is  by 
no  means  certain  that  the  applications  and  extensions  of  the  Kolmogorov 
hypothesis  to  turbulent  diffusion  are  firmly  based. 

The  main  purpose  of  the  present  Memorandum  is  to  examine  the  re¬ 
sults  of  applying  to  these  two  problems  --  the  diffusion  from  a  point 
source  and  the  spectrum  and  decay  of  homogeneous  fluctuations  --  a  trun¬ 
cated  Wiener-Hermite  expansion  of  the  velocity  and  concentration  fields. 

(9) 

This  method  has  recently  been  applied  by  Meecham  and  Jengv  to  the  de¬ 
cay  of  isotropic  turbulence  with  results  that  are  not  too  bad.  However, 
there  are  doubts  about  the  convergence  and  completeness  of  the  Wiener- 
Hermite  expansion,  and  it  has  recently  been  pointed  out^^ 


that  in 
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homogeneous  turbulence  the  truncation  prevents  energy  transfer  to  high 
wave  numbers  if  the  energy  spectrum  has  a  cutoff  at  the  initial  in¬ 
stant.  However,  the  present  writer  believes  that  the  objections  to 
the  method  arise  because  the  critics  take  literally  the  assertion  by 
the  methodfs  developers  that  turbulence  can  be  regarded  as  a  small 
perturbation  about  a  normally  distributed  velocity  field.  This  inter¬ 
pretation  is  clearly  untenable,  both  on  general  theoretical  grounds 
and  because  numerical  work  shows  that  terms  initially  small  do  not  re¬ 
main  small,  and  is  rightly  criticized.  On  the  other  hand,  the  method 
can  be  interpreted  in  the  spirit  of  the  integral  methods  used  in  ki¬ 
netic  theory,  boundary- layer  theory,  etc.  That  is,  a  finite  trunca¬ 
tion  of  the  Wiener-Hermite  expansion  can  be  regarded  as  an  approxima¬ 
tion  to  the  velocity  field,  the  approximation  being  determined  by  the 
requirement  that  it  satisfy  certain  integrals  of  the  Navier- Stokes 

equations,  in  a  manner  similar,  for  example,  to  the  Mott- Smith  or  Kar- 

•>v 

man-Pohlhausen  approximations.  It  goes  without  saying  that  the  re¬ 
sults  of  such  an  analysis  must  be  inconclusive,  since  in  the  complete 
absence  of  any  reliable  theory  there  is  nothing  with  which  they  can 
be  convincingly  compared.  But  it  is  precisely  because  the  prospects 
of  a  completely  sound  theory  are  so  distant  (and  the  experimental 
situation  is  far  from  satisfactory)  that  the  method  is  worth  pursuing. 
In  principle,  at  least,  the  accuracy  can  be  improved  by  taking  more 
terms;  a  virtue  that  some  of  the  other  crude  approximations  do  not 
possess.  Furthermore,  the  method  is  "honest"  in  the  sense  that  pre¬ 
conceived  physical  ideas  are  not  built  into  or  hidden  in  the  mathemat¬ 
ical  approximations. 


A  general  name  for  such  approximations,  which  we  henceforth  adopt, 
is  "Galerkin  methods." 
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TI.  THE  WIENER -HERMITE  EXPANSION 

The  basic  idea  of  the  method  is  due  to  Wiener,  and  was  ap- 

(12) 

plied  to  turbulence- type  problems  by  Meecham  and  Siegel^  and  others. 
The  technique  is  founded  on  the  white  noise  or  ideal  random  function 
a(x),  which  is  the  derivative  of  a  continuous  random  walk  and  has  the 
property  that  a(x)  is  a  normally  distributed  random  function  with  zero 
mean  and 


(a(x)a(x'))  =  6(x  -  x ') 


(4) 


One  then  considers  the  stationary  random  function  defined  by 


f ( x)  =  +  f  K^(x  -  xOh^^(x')  dx‘ 


f[  K^^(x  -  x  x  -  x/,)H^‘U(x/,  x"j  jjx  r  cix  + 


(5) 


where  the  Wiener -Hermite  functions  are  statistically  orthonormal  com¬ 
binations  of  a(x)  with 

=  1,  H^(x)  =  a(x),  H^(x,  x ')  =  a(x)a(x;)  -  6(x  -  xO,  etc. 

(6) 

The  kernel  functions  K  are  nonrandom  ordinary  functions.  It  has  been 

(13) 

shown  that  the  expansion  Eq.  (5)  is  complete,  but  the  question  of 

how  one  determines  in  practice  the  kernel  functions,  given  the  proba¬ 
bility  distribution  of  f,  seems  unanswered,  except  when  f  is  itself 
normally  distributed  and  the  problem  is  trivial.  If  f  satisfies  an 
equation  of  motion,  substitution  of  Eq*  (5)  and  equating  coefficients 
of  the  gives  equations  for  the  nonrandom  .  Since  in  general 

there  is  an  infinite  number  of  coupled  equations,  the  usual  difficulty 
of  turbulence  theory  is  not  avoided,  but  the  hope  is  that  approximate 
solutions  will  give  useful  results. 

For  the  problem  of  turbulent  diffusion,  we  shall  take  the  ideal 
random  function  to  be  a  random  vector  function  of  both  position  vector 
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x  and  time  t.  We  then  consider  the  normally  distributed  vector  func¬ 
tion  a_^(x,  t)  with  the  property 

^a^x,  t)a^(x,  t)  )  =  6i,^(x  -  x')6(t  -  t ')  (7) 

We  now  expand  the  random  concentration  9(x,  t)  as  a  series  of  statis¬ 
tically  orthonormal  Wiener-Hermi te  functions  obtained  by  generalizing 
Eq.  (6): 

9(x,  t)  =  ®(x,  t)  +  f  Ck(x,  t,  T)ak(§,  t)  d§dT 

+  j  Dkje(2E>  I>  I'>  T,  t') 

{ak(l,  t  ')  -  6kje6(i  -  £7)6(t  -  T  7)}d£d£7dTdT  7  +  ... 

(8) 

The  scalar,  vector,  and  tensor  kernels'®,  C^,  ,  etc.,  are  or¬ 

dinary  nonrandom  functions  of  their  arguments.  So  far  no  assumptions 
have  been  made  about  statistical  stationari ty .  Mean  values  should  be 
interpreted  as  ensemble  averages.  Clearly 

{9  (x,  t)^  =  ®(x,  t)  (9) 

However,  there  are  no  analogous  relations  for  mean  values  of  products 

of  0;  for  instance,  ^0^  involves  integrals  of  kernels  of  all  orders. 

In  this  respect  the  Wiener-Hermi te  expansion  is  clumsy.  Apparently 

it  has  not  been  proved  for  random  functions  of  more  than  one  variable 

(13) 

that  the  expansion  Eq.  (8)  exists  and  is  complete,  but  the  proof 

for  functions  of  one  variable  should  be  readily  generalizable . 

(9  12) 

The  main  applications  of  the  method  so  far  reported  ?  have 
not  used  the  expansion  Eq.  (8),  but  an  alternative  form  in  which  the 
ideal  random  function  is  independent  of  t  and  a  function  of  x  only. 

Now  for  the  problem  of  decaying  turbulence,  in  which  the  statistical 


-6- 


properties  are  given  at  t  =  0  and  one  tries  to  calculate  how  mean  val¬ 
ues  change  with  time,  it  would  appear  that  an  expansion  in  the  re¬ 
stricted  Wiener-Hermite  functionals  based  on  a^(x)  is  equally  appro¬ 
priate  and  mathematically  complete,  and  that  the  time  dependence  is 
adequately  described  by  the  variation  of  the  kernels  with  t.  The 
method  thus  suffers  in  this  case  from  a  lack  of  uniqueness.  Indeed, 
for  decaying  turbulence  where  the  time  dependence  is  imposed  by  the 
equations  of  motion  and  not  by  the  statistics  of  the  ensemble,  it 
would  appear  that  we  are  free  to  choose  ideal  random  functions  of  x 
with  an  arbitrary  dependence  on  t.  Thus  an  infinite  number  of  differ¬ 
ent  expansions  may  be  possible,  for  each  of  which  the  kernels  would 
have  a  different  t  dependence.  The  relationship  between  different  ex¬ 
pansions  has  not  been  examined  analytically.  It  is  reasonable  to  sup¬ 
pose  that  all  expansions  will  give  the  same  final  answers  if  the  in¬ 
finite  number  of  equations  is  solved  exactly.  However,  truncations 
of  the  expansions  will  give  different  results.  (Similar  remarks  can 
be  made  for  the  application  of  the  method  to  inhomogeneous  turbulent 
flows  which  are  statistically  stationary  in  time.)  Clearly,  we  would 
like  to  choose  an  expansion  for  which  the  Wiener-Hermite  expansion 
converges  as  rapidly  as  possible,  in  order  that  the  error  in  the  ; 
Galerkin  method  is  small.  It  is  not  known  how  such  a  choice  should 
be  made,  but  it  appears  that  choosing  a_^(x,  t)  to  be  a  stationary  ran¬ 
dom  function  of  t  is  less  restrictive  and  probably  "safer”  than  taking 
it  to  be  independent  of  t;  since  the  latter  choice  ignores  the  fact 
that  any  realization  will  drift  away  from  its  initial  state,  thus  re¬ 
quiring  more  and  more  higher-order  kernels  to  describe  it  in  terms  of 
the  initial  statistics  as  t  increases.  In  the  present  Memorandum  we 
shall  use  the  expansion  Eq.  (8)  with  ideal  random  functions  that  sat¬ 
isfy  Eq.  (7),  and  leave  for  further  study  the  question  of  whether  the 
results  are  better  or  worse  than  those  that  might  be  obtained  with  a 
different  time  dependence  of  the  ideal  random  function. 
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III.  DIFFUSION  FROM  A  POINT  SOURCE 


We  apply  the  method,  first  to  the  one-particle  diffusion  problem. 
We  examine  the  solutions  of  Eq.  (1)  with  u(x,  t)  a  given  stationary 
random  function  of  x  and  t,  and  the  initial  condition  that  at  t  =  0, 

9  =  ^00;  i-e.,  a  blob  is  released  at  the  origin  at  time  zero.  In  or¬ 
der  to  restrict  the  algebra,  we  shall  retain  in  the  Galerkin  method 
the  smallest  number  of  terms  that  gives  nontrivial  results.  Thus  we 
shall  assume  that  the  velocity  field  has  a  Gaussian  distribution  with 
zero  mean;  i . e . , 


u.(x3  t) 


K  (x  -  g,  t  -  T>a j  (§,  t)  dg  dT 


(10) 


where  K_  is  a  given  function.  For  the  concentration,  we  take  as  the 
approximation  to  i ts  distribution  the  first  two  terms, 


6(x,  t)  =  ®(x,  t)  +  !  C  (x,  g,  t,  T)a,(g,  t)  dg  dT 


(ID 


where  ©  and  are  to  be  found.  The  kernel  .  in  Eq.  (10)  is  a  func¬ 
tion  of  x  -  |  and  t  -  t  because  the  velocity  field  is  statistically 
stationary  in  both  x  and  t,  the  concentration  does  not  have  this  prop¬ 
erty  . 


Following  the  spirit  of  the  Galerkin  method,  we  substitute  Eqs . 
(10)  and  (11)  into  the  equation  of  motion,  Eq.  (1),  and  obtain  equa¬ 
tions  for  ©  and  C,  by  imposing  the  condition  that  the  result  of  the 

K. 

substitution  should  have  zero  mean  and  should  be  statistically  orthog¬ 
onal  to  a.(x>  t).  After  some  algebra,  we  obtain 


t)  + 


K.  .  (x 
ij  - 


5,  t  -  T)  ^  (x,  §,  t,  T)  dg  dT  =  K  ~  2  (x,  t) 


dx. 

i 


dx7 


(12) 


and 


-8- 


ac . 

(x,  §,  t,  t)  +  K  (x 


I,  t  -  t)  @(x,  t)  =  k  C  (x,  %,  t,  t) 

i  Sx.  J 

1 


The  initial  conditions  for  these  equations  are,  at  t  =  0, 


(13) 


0  =  6(x),  C  j  =  0 


(14) 


For  the  sake  of  simplicity  we  now  put  K  =  0,  and  thus  neglect  the 
effect  of  molecular  diffusivity.  For  diffusion  from  a  point  source, 
the  molecular  transport  is  generally  negligible.  Integrating  Eq.  (13) 
with  respect  to  t  and  differentiating  with  respect  to  x,  we  obtain 

SC. 

(x,  £,  t,  t) 
i 

Further, 

SK.  . 

3”1  (x  -  I,  t  -  t)  =  0  (16) 

i 

because  div  u  =  0.  Substitution  of  Eq .  (15)  into  Eq.  (12)  now  gives 
(with  h  =  0) 


!  Sx. 

O  1 


-  '*>, 


t)  t —  ®(x,  t  )  [  dt 

Sx.  —  j 
k 


05) 


ft® <5*  t} 


Sx.  , 

1  o 


f  dt '  |^-  (x,  t  0  Ki j  C2S  -  i>  t  -  T)Kkj(x  -  t'  -  t)  dj,  dT 


07) 


Now  From  Eq.  (10),  we  have 
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XT" 


R 


ik 


(t  -  t') 


=(  u_L(x,  t)uk(x,  t ') )  =  J  K_(x  -  5,  t  -  T)Kkj  (x  -  j>,  t  '  -  t)  d§  dx 

(18) 

where  is  the  Eulerian  one-point,  two-time  velocity  covariance 

tensor  (the  velocities  are  measured  in  a  frame  with  zero  mean  velocity). 
Substituting  Eq.  (18)  into  Eq.  (17),  we  obtain 


§ 


dx. 


idxk 


@(x,  t')Rik(t  -  t')  dt' 


(19) 


as  an  integro-dif f erential  equation  for  the  mean  concentration©. 

When  ©  is  found  from  Eq.  (19),  follows  from  Eq.  (15)  and  the 
structure  of  the  concentration  field  0  is  given.  Thus  the  Galerkin 
method  applied  to  the  Wiener-Hermi te  expansion  yields  relatively  sim¬ 
ple  equations  without  too  much  labor,  at  least  if  molecular  transport 
is  neglected.  The  question  remains  whether  the  results  are  reasonable. 
Let  T  be  the  time  scale  of  the  Eulerian  covariance.  It  could  be 

Ei 

defined,  for  instance,  by 

T  =  :  R,.(t)  dt/R  (0)  (20) 

Ei  <.  11  11 

O 


Then  for  t  »  T^  it  is  almost  obvious  that  Eq.  (19)  reduces  to 


where 


3© 

3t 


D. 


ik  3x.3x, 
l  k 


© 


D 


ik 


Vt} de 

O 


(21) 


(22) 
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Thus  for  large  time-values  the  diffusion  of  mean  concentration  is  given 
by  a  heat-diffusion  equation  whose  effective  diffusivity  is  given  by 


Eq.  (22).  This  result  is  consistent  with  that  obtained  by  Taylor 


(1) 


[see  Eq.  (3)]  as  generalized  to  three-dimensional  diffusion,  and  says 
that  the  integrated  Lagrangian  correlation  coefficient  is  given  by  the 
integrated  Eulerian  correlation  at  a  fixed  point  (i-e.,  at  a  point  mov¬ 
ing  with  the  mean  velocity).  This  conclusion  is  plausible  and  rein¬ 
forces  the  belief  that  the  Wiener-Hermi te  expansion  can  give  reasonable 
results . 

On  the  other  hand,  for  small  t  we  can  replace  -  t  /)  by 

Rik(0)  -  4.uk)  ,  and  on  differentiating  Eq.  (19)  we  obtain 


a2©  /  \ 


a2e 

dx.  5x. 
l  k 


(23) 


i.e.,  for  small  time-values  ©  satisfies  a  wave  equation.  For  the  sake 
of  illustration,  it  is  useful  to  consider  isotropic  turbulence  for 
which 


Rik(t>  ■  u  6ikRE(t> 


(24) 


Then  Eq.  (23)  becomes 


2  2^ 
u  V  © 


which  has  the  solution 


(25) 


©  =  ut)  (26) 

where  r  =  |x|  and  6  1  is  the  derivative  of  the  (one-dimensional)  Dirac 
6-function.  However,  the  solution  [Eq .  (26)]  is  physically  unreasonable, 

because  it  predicts  that  the  mean  concentration  lies  only  on  a  spheri¬ 
cal  surface  of  radius  ut,  and  also  gives  both  positive  and  negative 
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values.  Thus  for  small  time-values  the  Wiener-Hermi te  truncation  is 

* 

apparently  unreliable. 

It  is  worth  examining  this  difficulty  a  little  further  to  show 
that  it  is  not  a  consequence  of  the  approximation  that  led  to  the  wave 
equation  [Eq.  (23)].  We  take  Laplace  transforms  (denoted  by  a  tilde) 


®(x,  p)  e  pt®(x,  t)  dt 

o 


(27) 


etc.  Then  Eq.  (19)  becomes  (for  isotropic  turbulence) 


u2RE(p)V2©(x,  p)  -  p@(x,  p)  =  -6(x) 


(28) 


with  the  solution 


®(x,  p)  = 


1 


-kr 


4ttu2Re(p)  1 


k  = 


Lu2Re(P) 


1/2 


(29) 


As  P  -  00,  Re(p) 


-1 


,  and  it  follows  from  the  properties  of  La¬ 


place  transforms  that  ©  is  discontinuous  at  r  =  ut.  Thus  the  wave 

1/2  1/2 

property  of  Eq.  (19)  persists  for  all  time,  but  if  t  »  T  the 

1/2  h 

concentration  is  spread  out  over  a  distance  u(tT  )  ,  which  is  small 

E 

compared  with  the  cutoff  radius,  and  the  nonphysical  cutoff  is  not 
significant.  (But  strangely  enough,  casual  observation  shows  that 
the  edge  of  a  cloud  of  diffusing  material  is  sharp.) 

If  we  take,  for  the  sake  of  illustration. 


-t/T 

Vo  *  e  E-  ■  — — r  o°> 

p  +  te 


* 

This  difficulty  is  not  due  to  the  use  of  an  ideal  random  func¬ 
tion  that  fluctuates  in  time.  The  wave  equation  is  also  obtained  if 
one  uses  a  time-independent,  ideal  random  function  of  position. 
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we  find  after  some  analysis  that 


where  H  is  the  Heaviside  function  and  I  is  the  modified  Bessel  func- 

v 

tion  of  order  v.  This  expression  agrees  in  the  limits  t  -»  «>  and  t  -*  0 
with  the  values  obtained  by  approximating  the  equation,  and  confirms 
the  conclusions  reached  above. 

However,  the  inadequacy  of  the  solution  for  small  t  may  not  re¬ 
sult  from  a  general  failure  of  the  method,  but  rather  from  the  fact 
that  the  method  is  unsuitable  for  an  initial  concentration  in  which 
all  the  material  is  concentrated  at  one  point,  since  two  terms  of  the 
expansion  are  insufficient  to  describe  the  distribution  when  it  is 
highly  concentrated.  The  Galerkin  method  may  work  better  when  the 
material  is  initially  more  diffuse  or  is  closer  to  a  concentration  of 
a  type  that  exists  some  time  after  the  initial  conditions.  In  other 
words,  although  the  solution  of  Eq.  (19)  with  initial  condition  Eq. 
(14)  may  not  describe  adequately  for  small  t  the  real  diffusion  of  a 
point  source,  the  error  may  cancel  out  when  this  solution  is  used  as 
the  Green's  function  to  calculate  the  spread  of  an  initially  diffuse 
cloud . 

To  demonstrate  this  argument,  consider  the  exact  solution  for 
small  t : 


0(x,  t) 


l  5t 


+  0(t3) 


(32) 
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where  =  u(x,  0)  and  0(x,  0)  =  9^^  is  the  initial  concentration. 

Suppose  also  that  the  initial  concentration  is  statistically  uncorre¬ 
lated  with  the  velocity  field.  Expanding  Eq.  (32)  and  taking  mean  val¬ 
ues,  we  obtain 


®(x,  t) 


(“M0)) 


dx. dx . 
i-  3 


+  •  •  • 


=  ® 


(0) 


2  2 
u  t 


vV0) 


0(t3) 


(33) 


for  isotropic  turbulence.  Now  we  solve  the  problem  using  the  Wiener- 
Hermite  truncation;  i.e.,  we  solve  Eqs .  (12)  and  (13)  with  the  initial 
condition  8  =  .  For  small  time-values,  Eq.  (26)  gives  the  funda¬ 

mental  solution  (since  the  equations  for  9  are  linear),  and 


0(x,  t)  =  -  f  777  6  '(r  -  ut)©(§,  0)  d|  (34) 


where  x_  =  x  -  Expanding  3(j|,  0)  as  a  power  series  about  x  and  in¬ 

tegrating,  we  find  that 


®(x,  t) 


zl  f  A'C*  -  ut)  ;-@(0) 

4  n  J  r  L 


+  r  . 
1 


d© 


(0) 


1  d2@(0)  , 

2  rirj  9x.3x. 


dr 


=  © 


(0)  ,  1  2  2  2^(0) 
v  +  -  u  t  v  ' 


(35) 


Thus  for  small  time- values,  the  Wiener-Hermi te  expansion  and  the  exact 

2 

solution  agree  to  order  t  when  the  initial  distribution  is  continuous 
rather  than  concentrated  at  a  point. 

There  is  therefore  reason  to  believe  that  the  method  is  useful 
when  applied  with  care  and  in  appropriate  circumstances.  Since  it  is 
only  an  approximation,  it  is  unreasonable  to  expect  more,  or  that  it 
should  always  give  good  results  even  if  misused. 
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We  can  use  the  results  of  our  analysis  to  find  expressions  for 

/  /:  \ 

the  Lagrangian  auto-correlation  coefficient.  It  has  been  shown'* 
that 


ALct)  - 


\(  HjU,  t)u1(0,  0)9 (x,  t)y  dx 


(36) 


where  0  is  the  solution  of  the  diffusion  equation  with  0  =  6  (x)  at 
t  =  0;  for  simplicity  we  suppose  the  turbulence  is  isotropic.  Substi¬ 
tuting  the  expansion  Eqs.  (10)  and  (11)  into  Eq.  (36),  we  obtain 


u2RL(t)  =  [(  UjCx,  t)u1(0J  0) )  ®(x,  t)  dx  (37) 

where  ©  is  the  solution  of  Eq*  (19)  with  initial  condition  Eq.  (14). 
Thus,  given  the  two-point,  two-time  Eulerian  correlation,  Eq.  (37) 
yields  the  Lagrangian  auto-correlation.  For  consistency,  it  should 
be  possible  to  show  directly  from  Eq.  (37)  that 


.-,00 

I  Rt (t)  dt  = 

J  O 


00 

Rp(t)  dt 

Jo  E 


(38) 


but  this  problem  has  so  far  defied  attack- 

The  formula  [Eq .  (37)]  has  the  same  structure  as  a  hypothesis  some¬ 
times  called  "Corrsin’s  conjecture."^)  There  is  an  important  differ¬ 
ence,  however,  since  in  Corr sin's  conjecture  the  exact  mean  value 
should  be  substituted,  whereas  in  the  present  analysis  we  should  sub¬ 
stitute  the  approximate  value  of  0  given  by  the  analysis. 

For  the  problem  of  diffusion  from  a  point  source,  it  is  easy  to 
interpret  the  Wiener-Hermite  truncation  in  terms  of  approximations 
made  directly  to  the  convection  equation.  Thus,  if  we  write 


0(x,  t)  =  @(x,  t)  +  0(x,  t)  (39) 

Eqs.  (12)  and  (13)  are  obtained  on  substituting  the  complete,  nontrun- 
cated,  Wiener-Hermite  expansion  of  0  together  with  the  velocity  given 
by  Eq.  (10)  in  the  pair  of  equations 
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S® 

St 


d0 
\  Sx. 


2 

HV  @ 


(40) 


and 


30 

St 


30 

u*  ^ — 
x  3x. 
x 


2 

HV  0 


(41) 


Now  Eq.  (40)  is  exact,  and  is  obtained  by  taking  the  mean  of  Eq .  (1)  . 
Then  subtracting  Eq.  (40)  from  Eq.  (1)  gives 


30 

at 


u, 

X 


HV  0 


30 
Ui  3x 


30 

u.  — 

i  3x 


> 


(42) 


which  demonstrates  the  error  in  Eq.  (41).  In  other  words,  use  of  the 
Wiener-Hermi te  truncation  is  equivalent  to  neglecting  the  fluctuations 
in  the  convection  of  the  random  part  of  the  field,  and  essentially  re¬ 
places  u_^  30/3x^  by  its  mean  value.  For  h  =  0,  Eq.  (41)  can  be  inte¬ 
grated  with  respect  to  t,  and  substitution  into  Eq.  (40)  gives  Eq.  (19) 
'k 

directly.  It  is  curious  that  Eq.  (41),  which  seems  a  natural  approx¬ 
imation,  does  not  appear  to  have  been  used  before  in  the  present  con¬ 
text. 

✓V 

I  am  grateful  to  Dr.  S,  C.  Crow  for  pointing  out  this  alterna¬ 
tive,  and  indeed  more  general,  method  of  obtaining  Eq.  (19). 
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IV.  DECAY  OF  HOMOGENEOUS  FLUCTUATIONS 


We  now  apply  the  method  to  the  decay  of  statistically  homogeneous 
fluctuations  of  concentration  in  a  statistically  homogeneous  and  sta¬ 
tionary  velocity  field.  Again  we  retain  only  the  minimum  number  of 
terms  necessary  to  give  nontrivial  results.  The  velocity  field  is 
given  by  Eq.  (10).  For  the  concentration  field,  we  now  use 


0(x,  t)  =  [  C,(x  -  £,  t,  f)a,  (§,  t)  d| 


dT 


+  [  Dw,(2£  ■  1>  2  -  1>  T>  t>  t/)  dI  dHdTdT 


k  a 


(43) 


where 


n,  T,  t')  =  a  (i,  T)a  (H,  t')  -  6,  .6(1  -  Tp6(T  -  r ')  (44) 


Since  we  wish  mean  properties  to  be  independent  of  position,  we  suppose 
that  the  mean  concentration  is  zero  and  that  the  kernels  are  functions 
of  x  -  |  and  x  -  T[. 

We  substitute  Eqs .  (10)  and  (43)  into  the  diffusion  equation  and 
require  that  the  resulting  expression  be  statistically  orthogonal  to 
the  Wiener-Hermite  functionals  of  the  first  and  second  order.  After 
some  reduction,  and  use  of  the  relationship 


,(2) 


Hk£  (|>  t,  T')a.(6,  t")  = 


HkM  +  ak(-’  T)’‘ii50  ‘  *  T"> 


+  a^(T[,  "  T")  (4 5) 


we  find  that 
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3t  Ck(‘-  '  ts  t) 


-  H>  t  -  T  0  Dkj(i  ■  1»  J£  -  H»  t>  T>  T  0  d£  dT; 


K^Qs  ~  H>  t  -  T/)  Djk(x  -  H»  x  -  1>  t>  T  t)  dn  dr ' 


+  H  - -  Ck(x  -  t,  T) 

9x. 

i 


(46) 


Dkj^  -  i>  X  -  1>  C>  T,  T  0  =  -K  (x  -  t  -  t')  ~~  Ck(x  -  ?,  t,  t) 


+  H  — 2  DkiCi 

dx2  kj 

1 


-  5,  x  -  T| ,  t,  T,  T  ') 


(47) 


[There  is  some  arbitrariness  in  the  tensor  Dk^ ,  since,  as  is  clear 
from  Eq.  (43) ,  an  arbitrary  antisymmetrical  tensor  may  be  added  to  D,  . 

KJ 

without  changing  9.  It  appears  to  be  more  convenient  not  to  impose 

the  requirement  of  symmetry  and  to  define  D,  ,  by  Eq.  (47).] 

kj 

These  equations  for  the  unknown  kernel  functions  and  are 

not  simple.  To  evaluate  the  present  approach,  we  shall  assume  for  sim¬ 
plicity  that  h  =  0  and  that  the  variables  are  statistically  isotropic- 
Now,  by  change  of  variable  and  use  of  the  equation  of  continuity,  we 
can  rewrite  Eqs .  (46)  and  (47)  as 


^Ci(x,  t,  T) 


5 

dx, 


J  Kkj(I’  1  "  T^tDij^-’  C’  T’  T  ^  +  D j i (£>  c>  T  '>  T)}dS  dT  ‘ 


(48) 


respectively. 

Let  us  consider  now  the  distribution  at  time  t,  given  that  at  the 
initial  instant  t  =  0  the  non-Gaussian  part  of  the  concentration  was 
zero.  Then  Eq.  (49)  can  be  integrated  with  respect  to  t  and  the  result 
substituted  in  Eq.  (48).  After  some  reduction  we  obtain 

§£  Gi(x,.  t,  t)  =  V2  J  RE(t  -  t,)Ci(x,  t',  t)  dt /  (50) 

o 

where 

&i Ct  -  t  ')  =  (  u.(x,  t)u  (x,  t  ')  )  (51) 

is  the  Eulerian  one- point,  two-time  covariance,  and  we  have  used  Eq. 

(16)  and  the  principle  of  isotropy,  which  implies  that 

'Kkjd>  t  -  tO  g|~C.(_|,  t',  t  ')  d£,dr'  =  6kAF(t,  t') 

Notice  that  the  second  time- variable  t,  which  relates  the  concentration 
to  the  representations  of  the  ideal  random  functions  at  different  times, 
only  enters  parametrically  into  the  equation.  After  Eq.  (50)  has  been 
solved  for  C^,  Eq.  (49)  determines  the  D_ .  It  is  not  surprising  that 
the  Gaussian  part  of  the  concentration  field  satisfies  the  same  equa¬ 
tion  as  the  mean  concentration  in  the  diffusion  from  a  point-source 
problem.  The  remarks  made  in  Section  III  about  the  equation  still  ap¬ 
ply,  and  we  expect  the  results  of  the  present  analysis  to  have  some 
physical  significance. 

We  now  turn  to  the  correlation  and  spectrum  functions  of  9(x>  t). 
After  reduction,  it  is  found  from  Eq.  (43)  that 
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(  S(x)0(x  +  r))  =  S(r,  t) 


CjC?3  t>  T)Cj(j>  +  r,  t,  t)  d_|  dT 


+  [  JL>  t)  T’  T/^Djx^  +  -’  H  +  £’  c>  T>  T  Od^d^drdT  ' 


+  [  DU(I»  T)>  t,  T,  T')Df  i(l+  r,  §  +  r,  t,  t T )d^d]]_dTdT  * 

(52) 


jiN-’  ’  ’  '~Aj 

It  can  be  verified  from  Eqs .  (48)  and  (49)  that 


h s<0'  0  '  It  ^  ‘  0 


(53) 


Thus  the  truncated  equations  preserve  mean-square  concentrations  as 
they  should,  since  G  has  been  chosen  to  be  statistically  orthogonal 
to  D0/Dt. 

Expressing  the  right-hand  side  of  Eq*  (52)  in  terms  of  the  Fourier 
transforms,  we  find  that 


S(r,  t)  =  (2~)3  f  c.(k,  t,  T)C.(-k,  t,  r)e1-'r  dk  dr 
J  J  J 

+  (2tt)6  [  D  (k,  k',  t,  t,  tOD  C-k,  -k'»  t,  r,  T/)ei(-+-,),-Xdkdk/dTdT' 


+  (2tt) 6  f  D  (k,  k t,  T,  T')D  C-k',  -k,  t,  t  T)el(-+-/)  -dkdk'dTdx' 

(54) 

where 


C  (k,  t,  t  -  — — —  |  C.(x,  t,  r)e  *  dx 

J  (2tt)J  •’  J 


(55) 


5  ,(k,  k /,  t,  t,  t')  =  - -r 

(2n)6 


D^(x,  i>  t)  T.  T  )e - —  —  dxdj. 


(56) 
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The  equations  for  the  Fourier  transforms  are 

t,  t)  =  -ikp(2n)3  K^-k',  t  -  t  ') 

{^(k,  t>  T>  T  7)  +  ^(k',  k,  t,  t t)|  dk'dT,' 

(57) 

§7  ^(k,  -/j  t’  T’  T  ^  =  -ikp%i('-/’  t>  (58) 

from  which  we  deduce  that 

f^C  'Ck,  t,  t)  =  -k2  f  RE(t  -  t')c  (k,  t',  t)  dt'  (59) 

J  o  J 

If  we  define  a  spectrum  function  for  the  concentration  distribution  by 

A(k,  t)  =  — 1__  I  s(r,  t)e_it  L  dr  (60) 

(2tt) 

it  follows  from  Eq.  (54)  that 

4>(k,  t)  =  (2tt)3  C  (k,  t,  T)c\(-k,  t,  t)  dr 


+  (2tt)  6  f  D  (k  k  - 

J  j/~  - 

k t,  t, 

-k  +  k ,,  t, 

T,  t  ')  dk  'dTdT  ' 

+  (2")6 1  Vt'.  i  - 

1 ?r 

rr 

TOD^.C-k  + 

k',  -k',  t, 

t',  t)  dk/dTd-r' 

(61) 


For  an  isotropic  distribution  in  which  S  and  A  are  functions  of  the 
magnitudes  r  and  k  alone,  the  spectrum  function  with  respect  to  wave- 
number  magnitude  is 
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00 

r(k,  t)  =  4rrk^A(k,  t)  =  J*  s(r,  t)r  sin  kr  dr  (62) 


By  using  Laplace  transforms,  we  can  now  write  the  solution  of 

Eq.  (59)  in  closed  form.  It  is  easily  shown  that  C,  is  proportional 

pt  ,  ^ 

to  e  ,  where  p  is  the  root  of 


P  +  k  R£(p)  =  0 


(63) 


and  R  (p)  is  the  Laplace  transform  of  R  (t)  .  If  T  is  the  Eulerian 
L  EE 


time  scale,  then  it  is  clear  that 


C.(k,  t,  t)  =  C . (k,  0,  t) e"k  u  TEt}  for  k  «  (uT  )_1  (64) 
J  J  E 


and 


C  (k,  t,  t)  =  C  ^  (k ,  0,  t)  cos  kut,  for 


»  (uTe)_1  (65) 


The  wave  number  (uT  )  =  k  ,  say,  is  characteristic  of  the  ener- 

gy-containing  eddies  of  the  velocity  field.  This  wave  number  separates 
the  Gaussian  contribution  to  the  concentration  spectrum  into  two  parts- 
For  wave  numbers  less  than  k^,  the  interaction  causes  an  exponential 
decay  with  time;  for  k  »  k  ,  the  interaction  gives  a  rapidly  oscillat¬ 
ing  fine-scale  structure  to  the  spectrum.  The  oscillatory  behavior  is 
associated  with  the  wave- like  features  of  Eq-  (50),  and  may  be  a  spur¬ 
ious  effect  due  to  the  truncation.  However,  the  quantity  C always 
appears  squared  in  relevant  physical  quantities,  and  for  large  kut  it 

seems  appropriate  to  neglect  the  fine-scale  oscillatory  structure  and 
2 

replace  cos  kut  by  its  mean  value  of  one-half.  Then  if  T  (k,  t)  is 

G 

the  contribution  to  the  spectrum  function  from  the  Gaussian  part  of 
the  concentration  field, 
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2  2 

-2k  u  T  t 

r  (k,  t)  4=  r  (k,  0)e  for  k  «  k 


(66a) 


J  rG(k'  °> 


for  k  »  k  ,  kut  »  1  (66b) 

e 


The  exact  forms  depend  on  the  unspecified  Eulerian  correlation,  but  for 
the  sake  of  discussion  we  can  take  these  expressions  as  holding  respec¬ 
tively  for  k  <  k  and  k  >  k  . 

e  e 

We  consider  now  the  contributions  from  the  non-Gaussian  part  of 
the  concentration  field.  Quite  generally,  we  can  write 


Cj(k,  t,  t)  =  C. (k,  0,  T)C(k,  t) 


(67) 


Then  from  Eq.  (58), 


^  t 

D  (k,  k',  t,  T,  t  /)  =  -ik  C.(k,  0,  T)  K  ,  (k t'  -  t  /)C  (k,  t ')  dt  ' 
U  -  -  P  J  -  JQ  P£ 

(68) 

On  substituting  into  the  second  and  third  terms  of  Eq.  (54),  we  find 

after  some  reduction  (using  the  isotropy  and  incompressibility  of  the 

velocity  field)  a  contribution  S  (r,  t) ,  say; 

NG  — 


sng(£>  t) 


pt 


O  O 


k^k^A(k,  0)e^—  -^R^(£,  tr  -  t//)C(k;>  tOC(k,  t"  )dt  'dt^dk  (69) 


where 


R  (r,  t'  -  t")  =  (  u  (x,  t7)u  (x  +  r,  t")  \  (70) 

pq-  \  p  -  q-  -  / 

is  now  the  two-point,  two-time  Eulerian  velocity  covariance. 

The  expression  Eq.  (69)  can  be  given  a  simple  interpretation  if 
we  use  the  approximation 
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0 


£ 


o  o 


R  (r,  t'  -  t'^CCk,  t/)C(k,  t")  dt  'dt 
pq  - 


>a+»  = 


2k2u2  ><*  * 


R_  (r,  0)  (71) 


The  approximation  Eq.  (71)  is  in  fact  asymptotically  correct  for  large 

2 

k  if  we  use  Eq.  (65)  and  replace  cos  kut  by  one-half.  It  is  also  cor¬ 
rect  for  small  k  and  large  t  if  the  characteristic  time  scale  of 
Rpq(£>  t)  is  the  same  as  that  of  R^(t),  which  is  reasonable  for  small 
r.  Then  Eq.  (69)  can  be  written  approximately  as 


NG 


2u 


(r,  t)  =  — -  /u  (x,  t)u  (x+  r,  t) 


dr  dr  L 

p  q 


(r) 


(72) 


where 


V2  )  (r)  =  S(r,  0) 


(73) 


In  other  words,  the  non-Gaussian  contribution  to  the  concentration  co- 
variance  is  one-half  the  Eulerian  two-point,  one-time  velocity  correla¬ 
tion,  weighted  by  a  function  that  depends  on  the  initial  concentration 
covariance.  Note  that  the  Gaussian  contribution  is  given  by 


sG(r,  t) 


'  [C  (k,  t)]2A(k,  0)el-‘-  dk 


(74) 


and  it  follows  immediately  from  Eqs  .  (69),  (74)  and  (59)  that  S (0 ,  t) 
is  constant. 

Let  us  now  suppose  that  at  the  initial  instant  the  concentration 
spectrum  was  peaked  around  k^  and  was  negligibly  small  at  large  wave 
numbers.  Then  for  small  r  we  have  approximately 


S(r,  0) 


(75) 


where  X  k 
o  e 


1. 


It  now  follows  that  for  r  «  k 

e 


and  ut  »  r, 


-24- 
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S(r,  t) 


=  |  S(r ,  0)  + 


6 

pq 


r  r 


_  J£_a  _ 


r26 


10X 


?) 


We  can  rewrite  this  result  as 


(76) 


where  f(r)  is  the  longitudinal  velocity  correlation  function  of  iso¬ 
tropic  turbulence.  In  terms  of  the  concentration  spectrum  function 
F(k),  we  can  say  that  if  F  is  initially  Gaussian  and  peaked  on  the 
scale  of  the  energy-containing  eddies,  then  the  interactions  predicted 
by  a  Wiener-Hermi te  truncation  produce  at  larger  wave  numbers  a  spec¬ 
trum  due  entirely  to  the  non-Gaussian  contributions 

fi2 

r.(k)  =  — -  E(k)  for  k  »  k  (78) 

3uZ  6 


where  E(k)  is  the  energy  spectrum  function 

\  (F)  =  |  u2  =  J  E(k)  dk  (79) 

o 

-  2 

According  to  the  analysis,  the  transfer  of  spectral  density  of  9 
ceases  when  F  attains  the  value  given  by  Eq.  (78).  This  can  be  veri¬ 
fied  directly  by  calculating  the  transfer  func tion ^ 0 (x) 0 (x  +  ir)u^(x)  ^ 
and  showing  that  it  has  an  average  value  of  zero  as  kut  becomes  large. 
The  probability  distribution  does  not  attain  an  equilibrium  because 
the  Gaussian  contribution  continues  to  decay  with  a  slower  time  scale. 

If  the  spectrum  of  concentration  for  all  wave  numbers  were  exactly 
proportional  to  the  spectrum  of  energy,  we  would  have 
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?§2 

r(k)  =  ~2  E(k) 

3u 


Thus  the  analysis  has  predicted  that  for  k  »  k^  the  spectrum  of  con¬ 
centration  is  proportional  to  the  energy  spectrum,  but  is  only  half 
as  large  as  if  proportionality  existed  for  all  wave  numbers. 

In  the  next  section  we  shall  discuss  the  results  of  our  analysis 
and  compare  them  with  those  of  other  theories. 


-26- 


V. _ THE  SPECTRUM  OF  A  CONVECTED  SCALAR 


We  now  use  the  analytical  results  of  the  previous  section  to  build 
up  a  heuristic  theory  for  the  spectrum  of  a  convected  scalar  when  the 
Prandtl  number  is  not  small  compared  to  unity.  First,  it  is  necessary 
to  determine  the  reliability  of  the  Wiener-Hermite  truncation.  Al¬ 
though  we  do  not  expect  the  results  to  be  uniformly  valid  over  all 
length  scales,  it  is  reasonable  to  believe  that  they  may  describe  to 
a  good  approximation  some  of  the  properties  of  the  distribution.  We 
shall  therefore  make  the  hypothesis  (which  must  ultimately  be  tested 
by  experiment)  that  the  interaction  between  the  velocity  field  and  the 

concentration  field  is  well  described  by  the  analysis  for  scales  which 

3  1 ! 

are  larger  than  the  Kolmogorov  length  l  =  (v  /e)  ,  but  not  for  smal¬ 

ler  scales.  This  hypothesis  is  equivalent  to  saying  that  the  error 
introduced  by  the  truncation  is  not  important  for  wave  numbers  for 
which  the  energy  spectrum  is  not  exponentially  small,  and  therefore 
avoids  the  criticism  of  Grow  and  Canavan.  The  neglect  of  the  mo¬ 

lecular  diffusivity  h  is  also  justified  if  the  Prandtl  number  u/h  is 
not  small  compared  with  unity. 

PRANDTL  NUMBER  OF  ORDER  UNITY 

First  we  consider  the  case  where  the  Prandtl  number  is  of  order 
unity.  The  velocity  fluctuations  are  damped  out  on  a  scale  that  is 
smaller  than  l  by  the  viscous  stresses,  and  we  expect  similarly  that 
the  concentration  fluctuations  on  the  same  scale  are  likewise  damped 

by  the  conductivity.  Thus,  T(k)  should  decay  exponentially  for 

-1  -2 

k  >  l  •  The  rate  at  which  9  decays  is  given  by 

0° 

2 H  I  k2r(k)  dk  (80) 

J 


Then  using  Eq.  (78),  we  have  a  decay  rate 
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X 


2h9 


3u 


2  J 


k  E(k)  dk 


H§2e 

3uu2 


since 


e  =  2v  f  k2E(k)  dk 

J 

o 


(81) 


(82) 


The  idea  here  is  that  when  the  Prandtl  number  is  of  order  unity, 
concentration  fluctuations  are  destroyed  by  molecular  conduction  on 
the  scale  JL  •  These  fluctuations  are  produced  by  the  nonlinear  inter¬ 
action  of  the  velocity  and  concentration  fields  by  which  fluctuations 
with  a  scale  L  (that  of  the  energy- containing  eddies)  are  transformed 
into  a  smaller  scale.  We  assume  that  this  process  (which  can  be  called 
a  cascade,  although  we  do  not  imply  or  in  fact  believe  that  fluctua¬ 
tions  with  different  length  scales  are  statistically  independent)  is 
described  by  the  Wiener-Hermi te  truncation  analysis  with  molecular 
conduction  neglected,  according  to  which  the  spectrum  at  large  wave 
numbers  assumes  in  time  ^/u  (which  is  small)  the  form  (78);  and  that 
the  molecular  conduction  then  acts  to  reduce  the  values  of  this  spec¬ 
trum.  As  conduction  destroys  or  lowers  the  spectral  density,  the 
nonlinear  interactions  replace  the  spectral  density  to  retain  the  form 
(78)  for  the  spectrum,  where  9  is  now  a  function  of  time  and  has  a 
time  scale  large  compared  with  A/u.  In  other  words,  Eq  •  (78)  gives 
the  spectrum,  assuming  the  existence  of  a  quasi-equilibrium. 

Using  Eq.  (81),  we  can  write  Eq.  (78)  as 


r(k)  =  ~  E (k)  (83) 

0  H 

for  k  »  L  ^  and  v/k  of  order  unity. 

In  the  inertial  subrange  L  ^  «  k  «  i  \  the  unmodified  Kolmogorov 
theory  predicts  that 

-  Kc2/V5/3 


E(k) 


(84) 
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where  K  is  a  universal  constant.  Then  Eq.  (83)  gives 


r(k) 


.  aKs'1/3  V5/3 


(85) 


The  Obukhov-Corrsin  application  of  the  Kolmogorov  theory  (Ref.  4)  gives 


r(k)  =  K0xe" 1/3  k’5/3 


(86) 


The  result  of  Eq.  (85)  implies  that 

Kg  ^  K  (87) 

since  u/n  is  supposed  to  be  of  order  unity.  Thus  we  have  in  this  case 
confirmed  the  predictions  of  the  universal  equilibrium  theory. 

The  experiments  reported  in  the  literature^ ^ J were  carried 
out  for  Prandtl  numbers  large  compared  with  one,  and  cannot  be  used  to 
check  Eq .  (87) . 


PRANDTL  NUMBER  LARGE  COMPARED  WITH  UNITY 

(4) 

When  v /k  »  1,  we  can  apply  the  theory  of  Batchelor  to  derive 
the  spectrum  for  k  »  i  ^  or  the  correlation  function  for  r  «  X . 
Batchelor  showed  on  the  basis  of  a  physical  model  (which  makes  the  re¬ 
sult  far  more  convincing  than  if  the  argument  relied  solely  on  assump¬ 
tions  about  interactions  in  wave-number  space)  that  for  r  «  X,  the 
correlation  function  S(r,  t)  satisfies  the  equation 


X 


-yr 


dr 


2k 


s2s  i  ss’ \ 

.  2  r  Sr ) 

dr 


(88) 


where  y  (<  0)  is  an  effective  value  of  the  least  principal  rate  of 
strain  [see  also  Saffman^^].  Batchelor  assumed  that  the  straining 
field  responsible  for  the  fine-scale  structure  was  uniformly  distributed 
in  space  and  had  an  order  of  magnitude  given  by  the  root-mean-square 
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3f> 


1/2  1/2 

vortlcity  (e/v)  ;  a  plausible  value  of  y  is  then  -l/2(e/u)  •  How* 

ever,  there  are  strong  reasons  for  believing  that  the  vorticity  is  con¬ 
centrated  into  small  regions  in  which  its  order  of  magnitude  is  larger 
1/2 

than  (e/u)  .  The  small-scale  fluctuations  of  concentration  will  then 

also  be  concentrated  into  these  small  regions,  so  that  y  should  be  re¬ 
placed  by  a  quantity  somewhat  larger  than  in  Batchelor's  exposition. 

If  cr  is  the  intermi ttency  factor,  i.e.,  the  fraction  of  volume  in  which 
the  vorticity  fluctuations  are  large,  then 


(89) 


since  the  major  contributions  to  mean-square  vorticity  come  from  the 
regions  of  intense  shear.  Thus  it  is  appropriate  to  write  in  Eq.  (88) 


Y 


(90) 


The  value  of  a  is  open  to  question. 


Saf fman 


(8) 


has  argued  that 


a 


Re 


-1/2 


(V  /uL) 


1/2 


(91) 


and  in  the  following  discussion  we  shall  use  this  value;  but  it  should 
be  kept  in  mind  that  further  work  on  the  structure  of  turbulence  may 
suggest  a  different  value-  Equation  (88)  can  be  integrated^"*  and  gives, 
for  r  « 


S(r,  t)  = 


e2  .  * 

Y 


X  1  ,  1  yr2/4H 

-  — 1  r  e 
r  2 
o  ’  r 


e-yr  /4k  dr~j  dr  (92) 


-1  (4) 

The  corresponding  spectrum  function  for  k  »  X  is 


e*k  /y 

Yk 


r(k,  t) 


(93) 
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Thus  the  intermittency  which  increases  the  value  of  y  decreases  the  v 
amplitude  of  the  spectrum  function  in  the  "viscous-convective"  subrange 
and  also  decreases  the  size  of  the  conduction  cutoff 

je9  =  C-H/y)172  =  £(h/v)1/2  ct174  (94) 

Batchelor  pointed  out  that  Eq.  (93)  and  the  Obukhov-Corrsin  re¬ 
sult  (86)  are  continuous  at  k  =  JL  ^  if  a  and  are  both  of  order 
unity;  i.e.,  if  there  is  no  intermittency.  If  there  is  intermittency, 
there  would  seem  to  be  a  serious  conflict  with  the  Kolmogorov  hypothe- 
ses  as  applied  to  the  convection  problem. 

In  the  present  work,  we  suggest  that  the  intermittent  fine-scale 
structure  should  be  retained,  but  that  we  drop  the  Obukhov-Gorrsin  re¬ 
sult,  when  k  «  i  ,  replacing  it  by  the  results  of  the  truncated  Wiener- 
Hermite  expansion.  The  justification  is  similar  to  that  employed  above 
when  v /k  was  of  order  unity;  namely,  that  the  truncation  does  not  de¬ 
scribe  the  fine-scale  structure,  but  for  r  »  i  the  error  is  small  and 

we  can  regard  the  analysis  as  giving  a  quasi-steady  structure  for 

-  2 

L  »  r  »  the  slow  decrease  in  0  being  due  to  the  generation  of 
fine  structure.  Further,  the  small-scale  fluctuations  which  are  even¬ 
tually  dissipated  by  molecular  conduction  could  in  principle  be  analyzed 
by  retaining  many  terms  in  the  Wiener -Hermi te  expansion,  but  as  an  ap¬ 
proximation  are  given  by  Batchelor’s  theory  (modified  because  of  inter¬ 
mittency).  Thus  we  use  Eq.  (92)  or  Eq.  (93)  for  r  «  1  or  k  »  i  \ 
respectively,  in  the  so-called  viscous-convective  range;  and  Eq.  (77) 
or  (78)  for  r  »  1  or  k  «  l  ,  respectively,  in  the  so-called  inertial- 
convective  range. 

The  value  of  x  (which  is  not  assumed  here  a  priori  to  be  indepen¬ 
dent  of  the  Reynolds  number,  as  it  is  in  the  Obukhov-Corrsin  theory) 
can  be  found  by  matching  the  expressions  for  the  viscous-convective 
and  inertial- convective  ranges  at  r  =  1  or  k  =  i  However,  we  do 

run  into  a  troublesome  difficulty  which  has  not  been  resolved:  match¬ 
ing  in  physical  space  and  wave- number  space  gives  different  results. 

In  wave-number  space,  we  use  the  form  (84)  for  E(k),  since 
there  is  strong  experimental  support  for  this  expression.  For  Eqs . 


(93)  and  (78)  to  be  comparable  when  k  =  l  \  we  must  have 


82Ks 

X  “  3uV/2 


(95) 


It  follows  that  in  the  inertial  convective  range. 


rn  ,  1/2  -1/3  .  -5/3 

F(k)  sa  xCT  e  k 


(96) 


i.e. , 


K0  «  (Re) 


-1/4 


(97) 


if  we  use  Eq .  (91)  . 

These  results  are  now  different  from  those  of  the  Obukhov-Cor rsin 

theory.  Firstly,  x>  the  rate  of  decay  of  scalar  concentration,  is  de¬ 
l/4 

pendent  on  the  Reynolds  number,  behaving  like  (Re)  if  K  and  e  are 

independent  of  the  Reynolds  number.  Secondly,  the  coefficient  of 
-5/3 

Xk  decreases  with  increasing  Re,  and  is  not  constant  as  predicted 

by  the  arguments  based  on  the  hypotheses  of  the  universal  equilibrium 
theory.  Insofar  as  the  present  results  are  based  on  a  clear  mathemati¬ 
cal  analysis  rather  than  on  vague  assertions  about  the  statistical  in¬ 
dependence  of  Fourier  components,  they  deserve  to  be  taken  seriously. 

The  relatively  slow  dependence  of  Kft  on  Re  makes  it  hard  to  check  Eq. 
(97)  against  the  available  data,^^  especially  since  there  is  consider¬ 
able  uncertainty  about  the  evaluation  of  x* 

If  we  match  in  physical  space,  we  use  the  form  of  Eq .  (92)  for 
r  »  (-H/y)^^,  which  is^^ 


S  (r ,  t) 


(98) 


The  expression  for  r  »  1  is  in  physical  space 
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S(r,  t)  =  e2  -  1.2  Ks2/3  r2/3 

3u 


Equating  Eqs .  (98)  and  (99)  when  r  =  X  gives 


2.4  e2He  r.  ,  ,  1/2. --1 

X  ~  - 2  1/2  log  j 

3u  a  L 


(99) 


(100) 


This  does  not  agree  with  Eq.  (95),  the  difference  being  due  to  the 
logarithmic  term.  If  we  use  Eq.  (100),  we  would  have  in  the  inertial 
convective  range 

W1  N  1/2  -1/3  ,-5/3  -  ,  ,  1/2.  ,lni* 

r(k)  a  XCT  e  k  log  (v/hct  )  (101) 

and  the  spectrum  function  does  not  match  at  k  =  i 

The  difficulty  is  present  in  Batchelor’s  original  analysis,  and 
is  not  just  a  feature  of  the  current  approach.  Until  it  is  resolved, 
a  question  mark  hangs  over  the  assumption  that  we  can  determine  x  by 
matching  the  two  expressions  for  the  inertial-convective  and  viscous- 
convective  ranges  in  either  physical  or  wave- number  space.  In  any  real 
situation  the  difference  is  not  significant,  since  the  logarithmic  term 
is  not  large. 

A  crude  physical  justification  of  Eq.  (95)  can  be  given  to  demon¬ 
strate  that  it  is  consistent  with  intuitive  ideas  about  the  fine-scale 

structure.  We  believe  that  the  dissipation  is  due  to  fluctuations 

-2  1/2 

of  order  (0  )  having  length  scale  Xq  in  the  region  of  large  vortic- 
ity.  Then 


which  becomes,  with  the  use  of  Eq .  (94), 
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vQ2  1/2 
X  <=«  — r  a 

l 


and  on  using  the  hypothesis  that  5  =  vu  all  ,  we  arrive  at 


62e 

uV/2 


(102) 


Thus  the  Wiener- He mite  truncation  combined  with  Batchelor fs 
analysis  has  given  results  different  from  those  previously  obtained 
for  the  spectrum  of  a  convected  scalar  with  large  Prandtl  number.  It 
is  hoped  that  experimental  results  may  eventually  be  available  to  test 
the  theory.  In  any  event,  it  is  clear  that  progress  will  remain  dif¬ 
ficult  without  a  better  knowledge  of  the  fine-scale  structure  of  the 
velocity  field  in  turbulent  motion. 

In  this  last  connection,  it  is  perhaps  appropriate  to  mention 
here  that  the  modification  of  Batchelorfs  theory  to  take  account  of 
intermittency  does  not  affect  the  application  of  the  theory  to  the 
vorticity  field  of  turbulence  by  Saffman,^^  who  assumed  that  vorti- 

city  correlations  satisfied  an  equation  like  Eq .  (88) ,  with  Y  of 

1/2 

order  (e/v)  and  with  an  extra  term  to  take  account  of  vortex  line 

stretching.  If  there  is  intermittency ,  y  should  be  replaced  by  a 

value  characteristic  of  the  straining  field  distorting  the  vorticity. 

(8) 

For  reasons  explained  in  Saffman,  J  it  is  believed  that  this  straining 

field  is  not  that  of  the  vorticity  itself,  which  is  0(u/A),  but  rather 

that  of  the  straining  motion  whose  length  scale  is  the  Taylor  energy- 

2  1/2 

dissipation  parameter  A  =  (15v  u  /  e)  and  in  which  the  character- 

1  /2 

istic  rate  of  strain  is  u/A  ^  (e/v)  .  Thus,  changing  the  funda¬ 

mental  assumptions  does  not  in  this  case  change  the  formal  result. 
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